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ABSTRACT 

The experimental status of stable bound states made out of heavy quarks is reviewed. 
The need for a way to deal with the non-perturbative transitions involved calls for 
precision measurements on one hand, and for discovery of as yet undetected states to 
confirm predictions on the other. In this article, recent experimental contributions 
to heavy quarkonia spectroscopy and decay will be reviewed, mostly from data 
analyzed by the BES and CLEO collaborations. 

The most dramatic recent discoveries include the detection of the first 
stable L = 2 state, and a first non-pionic hadronic transition in the T system, and 
a first measurement of Xcj —* AA. Scans of the ip' and T(1,2,3S) resonances are 
used to add information on partial and total decay widths. 
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1 Heavy Quarkonia 



Heavy quarkonia are non-relativistic strongly bound quark-antiquark states. Ow- 
ing to asymptotic freedom, non-relativistic quantum mechanics should hold apply 
to a good approximation to heavy QQ systems [L. This allows us to compare 
heavy quarkonia to another bound particle-antiparticle state with slowly moving 
constituents, namely positronium. In fact, a similar spectrum of bound states is 
expected; see Figure^ The states are characterized by the following quantum num- 
bers: radial excitation n, total spin S, relative orbital momentum L, total spin 
J = L + S, parity P, and C . Not all states have been verified by experiment. Just 
as positronium, held together electromagnetically, provided a testing ground to gain 
understanding about that underlying force, studing heavy quarkonia allows us to 
take a better look at QCD. 

The large heavy quark masses result in small values of the running strong 
coupling constant in annihilation and production processes. In contrast to this, tran- 
sitions between QQ states, such as radiative de-excitation or splitting off a gluon 
pair that turns into a pion pair, thereby producing a lower-lying state, are soft pro- 
cesses. Thus many of the processes belong to the regime of non-perturbative QCD, 
a region that lacks thorough theoretical understanding as of now. 

The strong interaction has an impact on many measurements investigating 
weak physics as well. Also, new physics can turn out to be strongly coupled. This 
means that any opportunity that presents itself to study aspects of QCD as an 
example of a strongly coupled theory must be used. 

The bound ti state is very short lived since it can decay via the weak force, 
which is not possible for charmonium and bottomonium. This makes it unsuitable 
for studies of the strong interaction. Also, if the cc/bb states have masses above 
open their respective flavor production threshold (BB for bb states and DD states 
for cc systems), they decay fairly quickly. This leaves eight quasi-stable states in 
charmonium and thirty in bottomonium for spectroscopy studies. 

1.1 Theoretical Understanding 

As mentioned above, heavy onia spectroscopy highlights the soft regime of QCD, 
making it impossible to use perturbative methods to calculate transition quanti- 
ties. Perturbation theory does describe the long distance part of the heavy quark 
potential. 

This puts emphasis on the need for calculation methods that can handle 
non-perturbative phenomena such as Lattice QCD. Recent developments allow to 
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overcome the previous limitations that prevented unquenched calculations, i.e. that 
did not treat light quark loops, a neglect that was estimated to easily contribute 
uncertainties at the 10% level. As these become available, it is essential that their 
results be subject to critical comparison with experiment in order to enhance confi- 
dence before employing the methods to other areas such as the weak sector of heavy 
quark physics. 

Another solution is to use the quark velocity and a s as expansion param- 
eters, opening the way for effective theories of the strong interaction. 

Yet another approach is the use of phenomenological models, the parame- 
ters of which must be determined from experimental data. A common ansatz is a 
potential consisting of a Coulomb like term ~ 1/r, where r is the distance between 
the quarks, and a confinement term, linear in r. 
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Figure 1: The spectrum of stable bb states and allowed transitions within the system. 
Parallel arrows at the bottom stand for a symmetric configuation, producing a state 
of total spin S = 1, antiparallel arrows for an antisymmetric one (S = 0). Also 
indicated is the orbital momentum, increasing from left to right. A similar spectrum 
exists for charmonium. The states are labelled similarly, albeit with a subscript c 
instead ofb, and the 2S+1 Lj = 3 Si states are called ip (J/ip for n = 1). 
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Figure 2: Heavy quarkonia production diagrams. 

Production (left) and their corresponding decay (right) processes: 

a) e + e~ — > 7* — > QQ, b) 77 — > QQ; c j — > gluons — ► QQ- 

Quarkonium de-excitation by emission of two pions, e) creating charmonium from 
a B meson. 

2 Production 

Quarkonia can be produced in several ways, which reach different states within the 
spectrum. The first three listed here are mere reversals of QQ decay processes and 
are sketched in Figure b), and c). 

In electron positron colliders, the reaction e + e~ — > 7* —>■ qq results in 
states that can couple to a virtual photon, namely n 3 Si such as J/ip and T with a 
tiny admixture of n 3 D\. Direct resonance production offers the advantage of large 
production rates, giving access to branching fractions even as small as 10 -5 . 

Two-photon collisions allow direct creation of J = 0, 2 states, e.g. fj[ c ,b], 
X[c,6][o,2i- While they are readily available at e + e~ machines, they suffer from small 
production rates. Still they provide an important contribution in that they can be 
used for discovery purposes. 

Hadron machines, being able to produce any onia state in principle by 
annihilation of the pp pair into gluons, continue to contribute mostly to the study of 
production of charmonia. This environment suffers from large background; thereby 
one has to focus on exclusive decays. 

Two more scenarios: Downward transitions within the system provide an 
important route to otherwise not reachable states. Any collider that can produce 
B mesons, be it a hadron accelerator or an e + e _ machine running on the T(4S) res- 
onance, has access to cc states through weak decays of the b quark. These two 
processes are sketched in Figure Eld) and e). 

An important background for the reaction e + e~ — > QQ — > X, or more 
explicitly, e + e~ — > 7* QQ — > 7* — > X, is the case in which no intermediate 
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QQ resonance is formed. The presence of this channel adds to the measured cross- 
section both directly and by interference, which can be a sizeable contribution [2|. 
In most measurements, this contribution is not taken into account or subtracted. 
This background needs to be either measured, by running off the relevant resonance, 
or calculated. In measurements of the cross-section as function of energy (scans), 
the non-resonant production can be explicitly taken into account when fitting the 
line shape. 

3 Transitions 

3.1 Hadronic Transitions 

In order to conserve charge, transitions can either happen by emitting neutral parti- 
cles or a charged pion pair. Single n° transitions are isospin suppressed. Kaon pair 
transitions are phase space prohibited. 

3.1.1 Non-Dipion transitions 

By studying the decay ip' — > and plotting the invariant mass of the two 

photons, evidence for a single pion and 77 transitions can be obtained. The BES 
Collaboration, using their 15M ip' sample, will be able to perform precise measure- 
ments of these transitions. To illustrate the quality of the data sample, Figure |3] 
displays the result of a study of the decay ip' — > r y r yl + l~, where clear n° and rj signals 
are observed in the distribution of the invariant 77 mass [3]. 

\|r(2S) -> yyJA|/, J/v|/ ^ 1*1" y(2S) -> yyJA|/, JA|/ l + l" 

1000 F I J 
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Figure 3: Evidence for it and rj transitions in charmonium J3J/. 

The first non-pionic hadronic transition in the T system was recently found 
by the CLEO Collaboration. In their sample of (5.81 ±0.21) x 10 6 T(3S) decays evi- 
dence was found for the decay chain T(35) — > 7X&i,2 > jujT(1S) with lj — > ?t + /t~?t 
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and T(1S) — > l + l~ . (The decay x'm ~ * wT(lS) is phase space suppressed.) Re- 
quiring the presence of an T(IS') candidate, identified through its decay into a 
high-momentum lepton pair, guarantees that the background from udsc pair pro- 
duction in the data sample is negligible. The measured branching fractions, obtained 
by maximum likelihood fit to the energy spectrum of photons not assigned to the 
7T° candidates (Figure HJ, are at the percent level jl], thereby confirming the predic- 
tion that such transitions should have sizeable rates. 
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Figure 4: First observation of a non-pionic hadronic transition in bottomonium: 
Left: The three pion invariant mass, indicating the presence of an u. Right: Fit to 
the photon energy spectrum, with the individual contributions from T(3S) — > JXbi 
(dotted) and T(3S) — > 7x'&2 (dash-dotted) overlaid. 



3.1.2 Dipion transitions 

When comparing the dipion invariant mass spectrum in if>' — > ipiT + iT~, T(2S) — > 
T(15)7r+7r-, T(3S) T(2S)7r+7r-, and T(3S) -> T(lS)n + n-, the last has a dis- 
tinctive double-peak structure (Figure El right) [S]. Many models have been devel- 
oped that try and fit this behaviour, albeit lack of precision did not allow a clear 
discrimination between the models. CLEO has presented new preliminary data [7j. 
The structure is confirmed in both n + 7T~ and reactions, measured exclusively 
and inclusively, see Figure El 

4 Charmonium Singlet States 

The ground state of charmonium, n c , has been studied for some time. Theoretical 
predictions of the mass and width of its first radial excitation have so far been limited 
to potential model calculations and based on drawing an analogy of the 5 = pair 
(j]' c , n c ) to the S = 1 pair (ip', J/ip). An early measurement by Crystal Ball 
yielded a particularly low r]' c mass, thereby introducing a lower value for the mass 
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Figure 5: The dipion invariant mass spectrum in hadronic onia transitions. Left: 
new preliminary CLEO data, right: previous results from Crystal Ball 161 (ip' — > 
ipir~^Tc~~ ) and CLEO ]$J (T{2S) -> T(lS)7r+7r- ; T(35) -> T(2£)7r+7r- ; and 
T(35) -> T(15)7r+7r-;. ' 

splitting Am = m(rf c ) — m(r] c ). Recent measurements of CLEO and BaBar using 
two photon events and Belle studying B — > Kr)' c — > KK s K~n + decays jH] give ample 
indication towards a higher value, as does a Belle measurement of the spectrum of 
the system recoiling against the J/ip |5J in e + e~ — > J/ipX. At the same time, the 
rj c mass has been remeasured in these experiments. A precision measurement was 
performed by BES using radiative decays in a 58M J/ijj data set [TO] . A summary 
of experimental data on the r] c and 7^ masses can be found in Tabled Of particular 
interest to theorists will be the mass splitting between the two n — 1, L = states, 
m(r]' c ) — m^ip'), in comparison to the n = states, used to calculate the strength of 
the spin-spin interaction term in non-relativistic potential models. 

Many attempts have been made to measure the l x Pi state in charmonium, 
labeled h c [TH]. So far, the results are suggestive, but not conclusive. In jTE], the 
branching fraction B(ip' — > 7i°h c ) is predicted to be as large as 3.7 x 10~ 3 , with 
B(h c — > 7?7 C ) as large as 50%. Combined with a reasonable detection efficiency, even 
a moderate size t/>' data sample should be able to test this prediction. 
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Table 1: Experimental r] c and rf c mass measurements. 



experiment 


m(r] c ) (MeV) 


m(rf c ) (MeV) 


Crystal Ball 
Belle 2003 (exclusive) 
Belle 2003 (recoil) 
BaBar 2003 (77) 
CLEOII/II.V 2003 (77) 
BES 2003 {J/ip decay) 


2980 ± 8 H2] 
2979 ± 2{stat) 
2962 ± 13(stat) 
2983.3 ± 1.2(stat) ± 1.8(syst) 

2984.7 ±2.1 
2977.5 ± l.O(stat) ± 1.2(syst) 


3594 ± 5 H3| 
3654±6(stat) ±8(syst) 
3622 ± I2(stat) 
3632.2 ± 5.0(stat) ± l.S(syst) 
3642.7 ± 4.1 ±4.0(syst) 


Potential models Il4j 




3594 - 3629 



5 Search For New States in Bottomonium Using Radiative T(nS) decays 

In contrast to the situation in charmonium, no singlet bb state has been observed yet: 
the bottomonium ground state l 1 ^, r]b(lS), has not been seen (and neither have its 
excitations), which is also true for the h b states, n P\. Since they are not accessible 
in direct production, they are searched for by making use of photon transitions. 

CLEO studied magnetic dipole (Ml) transitions from the triplet 5 states 1 , 
using 4 x 10 6 T(3S*) and 3 x 10 6 T(2S) decays. The reactions examined in an inclusive 
photon spectrum analysis are: T(3S) — > r/^S)^, T(3S) — > ri^lS)^, T(2S) — > 
r] b (lS)j, and T(3S) — ► h b (lP)7i7i with h b — > 777^. A search window was defined 
based on predictions of the hyperfine splitting between the rj b (nS) and T(nS) states 
and the h b and Xb states. As a preliminary result [7] , no signal was found, but upper 
limits at the 10 -3 level in the corresponding search windows have been set on the 
branching fractions. This already rules out some of the models. Also, it states that 
a significant improvement in data sample size will be needed to establish a signal. 

5.1 L = 2 in Bottomonium 

The T(1D), a L = 2 state, is unique in that it is the only stable high-L state in heavy 
quarkonia as all others lie above open production threshold. Providing experimental 
evicendence is not only a matter of principle, but allows for discrimination amongst 
various theoretical models that were tuned on L = 0, 1 states. A second way in 
which such a state is special is that it decays preferentially electromagnetically, 
thereby is comparatively narrow. It was searched for by the CLEO Collaboration 

lr These allow AL = 0. Direct transitions have An = 0, whereas the ones with An > are 
called hindered. Hindered transitions, due to small wave function overlap made possible only by 
relativistic corrections, are suppressed, but benefit from a E 3 dependence in the transition rate, 
where E is the energy difference between the states. 
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in T(3S) decays in four-photon cascades: T(3S f ) — > 7X&(2P), Xb(2S) 7^(11?), 
T(l£>) — > 7%6 (IP), Xb(l-P) — ► t'^(I'S')- A signal of 6.8a significance is seen at 
10161.2 ± 1.6 MeV|17j. inconsistent with J = 3. Theory predicts dominance of the 
J = 2 state over J = 1 by a factor of about six |T8] . The new state is therefore 
assigned to be the Y(1 3 L> 2 ). 



Precise cross-section determinations around a resonance as function of the center- 
of-mass energy, often referred to as "scans", are an excellent tool to determine the 
resonance parameters with as little bias as possible. 

6.1 BES ip' scan 

The BES Collaboration has determined ip' resonance parameters m particular 
studying the reactions ip' — *> hadrons, ip' —* ii + ii^J/ip, and ip' — > fifi. The quantities 
determined in a simultaneous fit to these cross-sections were T tot , Y^, and Tj/^ )7r + 7! - , 
the derived quantities are r^, B had , B h and B J / 1 p n+n -: 

• The channel if)' — > fifi can be combined with other leptonic width measure- 
ments, testing the sequential lepton hypothesis, which states that : 



should be the same number for I = e, /i, r. The denominator is about unity for 
I — e,fi and 0.39 for I = r. This BES determination of the muonic branching 
fraction yielded B^ = (9.2 ± 0.8) x 10~ 3 . The E760 Collaboration published 
a value of B e = (8.3 ± 0.5(stat) ± 0.7 (syst)) x 10 -3 . A separate (and first 
direct) measurement of B T was performed as well by the BES collaboration [20J, 
resulting in a value of B T = (2.71 ± 0.43 ± 0.55) x 10~ 3 . After applying the 
above correction factor to the the r partial decay width, it can be compared 
with the other leptonic branching fractions. Agreement within experimental 
errors is observed. 

• The ijj' — > tt + 7t~ J '/ip line shape is of importance since this decay is frequently 
used as a normalizing mode. The precision obtained on the branching fraction 
in this measurement is 4.4%, improved over the PDG2002 value of 5.2%. 

• Finally, using the relations Y tot = T had + T M + Y e + Y T and Y e = Y ^ = 
r r /0.39 and measuring the ip' hadrons cross section, Y tot was obtained as 



6 Scans 
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Tf ot = (264 ± 27)keV. This can be compared with a PDG2002 value of vf ot = 
(300±25)keV [TT] and app scan result by E760 of Tf ot = (306±36±16) keV EH- 

6.2 CLEO T(l, 2, 3S) scans 

CLEO hopes to improve the precision of the leptonic width T e from currently 
2, 4, 9% for T(1,2,3S') down to the level of 2% for each of the three. One reason 
is that this parameter enters many other measurements. In addition, it provides a 
high-precision test for Lattice QCD, which begins to be able to reach this level of ac- 
curacy. Preliminary results show a statistical precision of 0.1, 0.3, 0.5%. Systematic 
errors are still being evaluated. 

7 Decays 

Heavy onia can decay via the electromagnetic or the strong force. Possible decay 
mechanisms for a heavy onium state are annihilation of the two heavy quarks into 
two leptons, one or more photons, or two or three gluons. For states below open 
flavour production threshold, the decay into leptons only contributes little to the 
total width (12% for the J/ip, less for other states), and the remaining hadronic rate 
is by far not accounted for by the exlusive decays measured so far. To some, the only 
pieces of information of interest are the leptonic decay widths, since the rest has to 
consist of radiative or hadronic decays. How exactly this greater remainder of the 
total rate is divided up into specific final states is largely governed by fragmentation 
dynamics. Electromagnetic decays have been discussed in Section I6~T1 

7.1 First Evidence of Xcj — > A A 

BES has reported the observation of Xcj ~^ AA Besides this being the first 

measurement of the branching fraction, this decay is of interest in comparison to 
Xcj PP- It has been shown that the lowest Fock-state expansion of charmonium 
states ( "color singlet model" , "CSM" ) is insufficient to describe P-wave quarkonium 
decays, both inclusively and exclusively, and that use of the next higher Fock-state 
("color octet mechanism", "COM") improves the agreement with experiment. The 
agreement of COM based prediction with the total measured width of the Xco as 
well as that for the partial width of Xcj — ► PP, obtained by using a carefully tuned 
nucleon wave function 22J, was encouraging. Generalizing the nucleon wave function 
to other baryons lead to a prediction for the partial width of Xcj ~^ AA as being 
about half of that for Xcj — ► PP |22j for J = 1,2. 
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Table 2: Experimental results on Xcj ~^ AA ; and comparison with Xcj —* PP- 





B( Xc j - AA) in 10" 4 


B(Xcj -> PP) in 10" 4 


J = 

J = l 

J = 2 


4.7ii;|±1.0 
2.6±5;g ± 0.6 
3.3±1;| ± 0.7 


2.2 ±0.5 
0.7 ±0.3 
0.7 ±0.1 



The analysis uses a recently collected 15M sample of if)' events, decaying 
radiatively to the Xc states. The desired channel is identified through one ir~p and 
one charge conjugated candidate. The distribution of the resulting AA events is 
shown in Figure El overlaid with the fit result, with the masses of the three Xc states 
as fit parameters. As a normalizing channel, ip' — > imJ/ip is used. The branching 
fractions thus derived are listed in Table |2l together with the result for Xcj —* PP 
from [IT]. The large experimental errors do not justify ruling out the prediction. 
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Figure 6: Distribution of the invariant AA mass in ip' — > 7(^77 )(p7r + ) events 123)/ . 



7.2 Hadronic Decays 

Predictions exist that relate the branching fraction for the l 3 Si into hadronic final 
states to those of its first radial excitation: if this decay happens predominantly 
via annihilation into gluons, then the wave function overlap should be the only only 
difference between the decay of the two states (aside from the difference in center-of- 
mass energy that the gluons have, which is not vastly different). This in turn can be 
taken from the leptonic branching fraction, another annihilation process. Ignoring 
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the running of the strong coupling constant 2 a s , one obtains 

_ B(i/>(2S) - H) ^ Bty(2S) - e+e~) 
V B(J/V> — H) B(J/ip e+e~) ' 1 J 

Using the leptonic branching fractions B(J/ip — > e + e~) = (5.93 ±0.1) x 10 -2 and 
B(ip(2S) -> e+e - ) = (7.3 ± 0.4) x 10~ 3 |TT] , the expected value for the ratio is 
Q = (12.3 ± 0.9)%. Different views exist as to whether this prediction is valid only 
for the inclusive process ip(nS) -*Ior also for specific final states. Some theorists 
posit that it should hold also for radiative decays. Moreover, a similar relation 
should exist for the ratio for the S = states. 

A number of channels have been studied in charmonium, including decays 
into a vector (V) and a pseudoscalar (P) particle, axialvector pseudoscalar, vector 
plus a tensor (T), radiative decays, multibody decays, as well as dibaryon final 
states |2H] • While many of these are not in outrageous disagreement with the 
above prediction, two VP modes are well-known for failing: prr (Q p7T < 0.65%) and 
K*K (Qk*k < 1%)- Also in VT states substantially lower ratios have been observed 
(Qujf 2 < 3.5%, Q pa , 2 < 2%). The overall picture is not clear, in part due to inaccurate 
experimental results. Many a theorist has given input on this anomaly that is 
often referred to as the pre puzzle, but so far none is able to explain all experimental 
results. It is of major interest to decide whether the J/ip rate is enhanced or the 
rate is suppressed. 

In the T system, a similar relation is expected to hold. As in this case two 
excitations are below dissociation threshold, two such ratios can be built. Using the 
corresponding leptonic branching ratios, one obtains 48% for T(2S) : T(1S) and 72% 
for T(3S) : T(1S). It is by no means clear what absolute rate to expect for T decays 
when extrapolating from charmonium. Depending on the model chosen to explain 
the p7r anomaly, the rates vary considerably. A preliminary CLEO result studying 
a variety of two-body hadronic T decays are upper limits of 4 • 10~ 5 or better on p7r, 
K*(892)K, pa 2 (1320), K* (892) K%( 1430), w/ 2 (1270), &i(1235)tt, and Ki{U00)K on 
the T(l,2,35') resonances ^7\. In particular, it is found that B(T(1S) — > pn) is 
well below 10~ 5 . 

8 Summary 

Heavy quarkonia continue to provide a testing ground for QCD calculations. 

A wide variety of measurements is being carried out, with hopes for many 
more results once BES and CLEO have analyzed their recently taken large datasets. 

2 A value of 0.85 is quoted for (a s (m^ 2 s))/a s (mj^)) 3 in 24 . 
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The goals for the future are to establish the states that have been predicted 
to exist but not observed yet and to provide precision measurement to allow a 
detailed comparison with theory. 

This will hopefully make it possible to gain further insight into the non- 
perturbative realm of QCD, from which many measurements in other areas will 
benefit greatly. 
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